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Abstract 

The supersymmetric formulation of t-J model is studied in this paper at 

the mean-held level where the 5-T phase diagram is computed. We find 

that slave-fermion-like spiral phase is stable at low doping concentration, and 

the slave-boson-like d-wave fermionic spin pairing state becomes energetically 

favourable when 5 > 0.23. An improvement in free energy using Gutzwiller's 

method lowers the transition doping concentration to 5 ~ 0.06. We also point 

out the existence of new branches of excitations in the supersymmetric theory. 
PACS numbers: 71.10.Fd, 74.25.Dw, 75.10.Jm 
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A lot of interest in studying the t-J model is raised since the discovery of cuprate 



mean-field theories which treat the constraint of no double occupany on average. The 
most successful mean-field approaches to the t-J model seem to be based on either the 
slave-fermion mean-field theory (SFMFT)i, which is successful at very small doping when 
antiferromagnetic correlation is important, or the slave-boson mean-field theory (SBMFT)i, 
which is successful at larger value of doping when the system becomes superconducting. The 
only difference between the two approaches is that two different representations of spin and 
hole operators are used in the two theories. Recently the focus has turned to the underdoped 
regime where crossover from antiferromagnetism to superconductivity takes placei'i. To un- 
derstand the complicated behaviour in this regime, it seems that a unified approach which 
incorporates the advantages of both SFMFT and SBMFT is essential. In a recent paper!, we 
show that it is in general possible to formulate lattice models with constraints of no double 
occupancy in such a way that the advantages of the two slave-particle representations are 
retained. The essential idea of the approach is to consider an enlarged Hilbert space where 
both slave-boson and slave-fermion representation of spins and holes coexist as independent 
states at each lattice site. A Hamiltonian is chosen such that the system in the enlarged 
Hilbert space is equivalent to 2^ {N = number of lattice sites) replicas of the original lattice 
model (supersymmetric formulation)!. 

In this paper, we present our mean-field calculation results on the supersymmetric for- 
mulation of t-J model. First we give some mathematical details of the theory. A spin-a 
state is represented by either 6^|0) or /J^|0), whereas a hole is represented by either JF"I"|0) 
or i3^|0), where are bosonic creation operators and f^, T'^ are fermionic ones. The 

Hamiltonian which is equivalent to 2^ replicas of t-J model in the Hilbert space spanned 
by the above states is! 



superconductorso. So far, analytical understandings of the model 




were largely based on 



E = -tY. {-hlj^,J^]h,^ + 6L + flB^J'lh,, + flB,B]f,, + h.c.) 
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The constraint of no double occupancy /ij/ij + Z^o-V'jWio- = 1 is imposed by a Lagrange 
multipher term as usual. The Heisenberg interaction can be re-written as 
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Af/ = ai^(3ji - anpj^, xtf 



J^aCilaPj^y where aio-,Pia can be either bi„ or fi„. Notice 



that only /^^^^^'^ and xTj'^ " terms appear in usual slave-fermion(boson) mean-field theory, 
whereas Aj^ and Xij ^^re supermatrices in our formulation. It is straightforward to show from 
Eq. that Hj is invariant under the super-unitary transformation hi Uihi^ipu —>■ Uiipa, 
where Ui is a local 2x2 unitary super-matrix. The hopping term can also be re-written as 



bbiff) 



Ht = -tY^ strinljXij + h.c. 

(id) 
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where 



n 



V 
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Notice that Ht is not invariant under the super-unitary transformation Ui because of the 
sign difference in the term —bl^J^iJ^jbj^ in Eq.(|I|). The sign difference is well known in 
SFMFtI. Thus although the Heisenberg model is supersymmetric, the t-term breaks 



supersymmetry. The Lagrangian L of t-J model in the enlarged Hilbert space is £ = 

Our mean-field theory is obtained by replacing Aj — * (Aj) = A, and by decoupling the 
quartic terms appearing in Ht and Hj. After some straightforward algebra, we obtain 

-t J2 str(7i:*,Xii + %Xij - T-CijXij + cc.) - ^ E str(x*jXii + XijXij - XijXij) 
-J E strA(A*,A,, + A^A., - A^A,,) (8) 

where Ajj(A*j), Xijixtj) ^'^d are Hubbard-Stratonvich type super- matrix field in- 

troduced to decouple the quartic terms in H. 

Next we integrate out hole and spin fields h and ip^^, arriving in an action in terms of 
the auxiliary super-matrix fields only. The supersymmetric mean-field theory is obtained 
by replacing the auxiliary super-matrix fields by their mean values which are determined by 
the saddle-point conditions, 

6S 6S 6S 



SQ^J{T) 6QUt) 5A,(r) 
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where Qij can be A°^, x'^f or Hif. Since the first-derivative of an action S with respect 
to any Grassmann field is a Grassmann variable and mean-field value of any Grassmann 
variable is zero, all the off-diagonal elements of the super-matrices vanish at the mean-field 
level and only diagonal terms of the supermatrix fields survive. Notice that both slave-boson 
and slave-fermion like mean-field parameters are present in the mean-field Hamiltonian. The 
relative weight of the two kinds of terms are determined by minimizing the free energy of the 
system with the constraint that the correct average number of spins and holes are obtained^. 

In our calculation we have only searched for mean-field solutions which preserve trans- 
lational invariance and time-reversal symmetry. Notice that the ground state energy of the 
Heisenberg term obtained from the corresponding variational wavefunction methodBEl is 
^var ^ 2EJ^, where Ej^ is the energy in our mean-field theory. This discrepancy is due to 



the defect of Hubbard- Stratonovich method when decouphng in more than one channelQ. To 
compare our results with those obtained from variational wavefunction method, the coupling 
constant J in our mean-field theory is scaled to 2J in the following. 

Fig.l shows the 5-T phase diagram with t/J = 3. The phase diagram can be divided 
into four regions. There are three low temperature phases at different values of 5 and a 
high temperature phase (IV) where all order parameters vanish. At low doping the stable 
low temperature phase is slave-fermion like (I), with the system characterized by antifer- 
romagnetic or spiral order, whereas at higher doping the stable low temperature phase is 
slave-boson like, with the system characterized by d-wave RVB pairing (II) or is a Fermi 
liquid (III) (at high doping). The phase transitions are all first-order (solid lines) except that 
between d-wave fermionic RVB and Fermi liquid (dot-dashed line) which is second-order. 
Notice that first-order transition to high-temperature phase is also obtained in pure slave- 
boson or slave-fermion treatment, and is believed to be a defect of mean-field theory. The 
same is also expected here. It is however, not so clear whether the first-order transition from 
slave-fermion-like to slave-boson-like solution is a defect of mean-field theory, or whether it 
really exists in J model. It is interesting to point out the similarity of our mean-field phase 
diagram to that obtained from SO (5) theoryH. Assuming that the first-order phase transi- 
tion to high-temperature phases are replaced by smooth crossover, the two phase diagram 
are qualitatively similar, except the presence of spin-charge separation and the replacement 
of antiferromagnetic phase by spiral phase in our mean-field theory. Whether there exists 
a hidden S0(5) symmetry in our theory is not clear at present. In the following, we shall 
describe the structure of the phase diagram in more details. 

At half-filling, antiferromagnetic state with only (A^^^) 7^ is the most stable mean-field 
solution until at sufficient high temperature T > 0.6 J, where (A^^^) vanishes and the system is 
in the high temperature phase similar to what happen in the usual SFMFT. For < 5 < 0.23, 
(Xij) becomes non-zero at low temperature, corresponding to the development of spiral order. 
Upon further doping the spiral phase becomes energetically unfavorable and is replaced via 
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a first-order transition by the slave-boson-like d-wave RVB pairing phase where only (Ajj) 
and {Xij) are nonzero. At still larger doping or higher temperature, the slave-boson-like 
d-wave RVB pairing phase is replaced by the Fermi liquid and "strange metal" phases, 
respectively with only (xf/) 7^ 0. Notice also that as in the usual SBMFT, a nonzero d-wave 
RVB pairing order parameter (AfJ) may describe either superconducting or spin-gap phases, 
depending on whether the holons are bose-condensed. The bose-condensation temperature 
Tb is found to be very high in supersymmetric mean-field theory as in the usual SBMFT 
(dashed line). Tb is believed to be suppressed by gauge fiuctuations once we go beyond 
mean-field theorylli. We have also searched for "mixed" states where both slave- fermion 
and slave-boson type order parameters coexist. However we found numerically that such 
phases are never the most stable. The most favorable phase is always characterised by slave- 
fermion-like or slave-boson-like order parameters except at high enough temperature where 
all the order parameters vanish. 

Notice that although the mean-field solutions we obtained are always either slave-boson- 
like or slave- fermion-like, there is an important difference between supersymmetric theory 
and usual slave-fermion or slave-boson approaches. In usual slave-fermion or slave-boson 
mean-field theories the statistics of the spin and hole excitations are fixed and are deter- 
mined by the types of mean-field theory being used. However, in the supersymmetric theory 
spin and hole excitations with both fermi- and bose- statistics existil For example, in the 
spiral phase where the ground state is the same as that obtained in the usual slave-fermion 
mean-field theory, we find that besides the usual bosonic spin-wave excitation, new disper- 
sionless high energy fermionic spin excitation of energy ~ 2.3J and bosonic hole excitation 
of energy ~ 106t appears in the excitation spectrum. The fermionic spin excitation energy 
decreases as the antiferromagnetic correlation weakens either because of raising temperature 
or increasing doping concentration, until the slave-boson-like d-wave RVB pairing phase be- 
comes energetically more stable. The excitation spectrum in the slave-boson-like phase has 
a "dual" structure to the antiferromagnetic phase. The low energy spin excitations are 
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fermionic as in the usual SBMFT. However, a dispersionless high energy bosonic spin ex- 
citation of energy ~ 0.4 J and fermionic hole excitation of energy ~ 1.5t (at low doping 
regime) are found to exist in the supersymmetric theory. Notice that finite dispersion is 
expected to appear in these high energy excitations after the fiuctuation arising from the 
off-diagonal elements of super-matrices is considered. 

One of us proposed that these high energy spin and charge excitations are in fact topologi- 
cal excitations in SFMFT and SBMFT, respectivelySi. It is suggested that the slave-fermion 
(boson) mean-field state can be considered as a topologically-disordered state of slave-boson 
(fermion) mean-field theory. In the transition from one mean-field state to another, the (low 
energy) elementary excitations of one state become the (high energy) topological excitations 
of another. The supersymmetric mean-field theory seems to support this picture. 

Notice that because free energy is an increasing function of 5 when 6 is small in SFMFli, 
there exists an instability region where the system is unstable towards phase separation (the 
area under dotted line) in mean-field theory. A possible consequence of this instability is 
that stripe phase might develop in this part of the phase diagram. Notice however, that 
our mean-field free energy is probably not accurate enough to describe the instability region 
correctly. 

The first-order phase transition between slave-boson-like and slave-fermion-like phases 
implies that our phase diagram is extremely sensitive to the free energy of the two phases. 
In particular, it is well known that the mean-field free energies are unreliable because of 
the relaxation of constraint of no double occupancy. A better way to compute free energies 
is to use Gutzwiller's projected wavefunction0 = P\(f)o), where |0o) is a mean- field 
state in the supersymmetric theory, and P = Ili6{hfhi + J^a 4'ia4'i(T ~ 1) is the projection 
operator. We shall estimate the energies of the wavefunctions {(pc) using an approximate 
method originated by Gutzwillei0 in the following. 

Let A be any operator where we want to find its ground state expectation value 
I^I^g)- Gutzwiller introduces the factor such thal0 (0g|A|0g) = fi'A(0o|^|0o) where 



the coefficient qa is determined approximately by considering only the ratio of classical 



weights of the Gutzwiller projected state and that of mean-field statellHO. Applying the 
above procedure to our supersymmetric mean-field theory, we find that there are four kinds 
of g factors, where two come from Heisenberg interaction term in slave-boson and slave- 
fermion mean-field states, respectively and the other two come from the kinetic energy 
terms. We found that gf = 4/(2 — for the Heisenberg interaction term in slave-boson 
mean-field state, where 6f = J^aifLfitr) = 1 in the slave-fermion mean-field state. 

Similarly, 

, ^ 2(1 + Sb) f ^ 2 + S, 

' (2 - v^ii^rr 2(1 - s^){En=i v^i^,ry 

for the hopping terms, where Sb = J^ai^la^ia) , Sjr = (J-'jj-'i) and Sq = {B\Bi). The details of 
these calculations will be reported in a separate paper. 

After applying the Gutzwiller's method, the resulting S-T phase diagram for the same 
value of t/J=3 is shown in Fig.2. Comparing with the bare mean-field phase diagram we 
find that the slave-fermion- like spiral phase has shrinked into lower doping regime. The 
transition doping concentration at zero temperature drops from 0.23 to 0.06, which shows 
that the importance of constraint in obtaining the correct free energies. The d-wave RVB 
pairing phase also extends to higher temperature region and covers a larger area in the phase 
diagram, in agreement with general experimental findings. 

In summary, we present in this paper the mean-field phase diagram of t-J model in 
supersymmetric representation. The supersymmetric theory produces several new features 
not present in usual slave-boson or slave-fermion treatment of t-J model. Besides providing 
a better phase diagram compared with experiment, new branches of spin and hole excita- 
tions are obtained in supersymmetric theory. The transition from slave-fermion-like spiral 
state to slave-boson-like d-wave superconducting state is first-order in our mean-field treat- 
ment, similar to what is obtained in S0(5) theory. The stability of our mean-field phase 
diagram against quantum fluctuations is unknown at present, and shall be a subject of our 
future works. Nevertheless, our success in constructing a theory which incorporates both 
the advantages of the slave-fermion and slave-boson mean-field theories suggests that super- 
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symmetric formulation can be a fruitful starting point to understand the relation between 
antiferromagnetism and superconductivity in t-J model. 
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FIGURES 

FIG. 1. Phase Diagram of J model with t/J=3. The phase diagram has four regions. Region 
I: slavc-fcrmion-like spiral phase. Region II: Slavc-boson-like d-wavc RVB pairing phase. Region 
III: Slave-boson-like normal state. Region IV: High-temperature phase. The area under the dotted 
line is unstable to phase separation in mean-field theory. The dashed line is the bose-condensation 
temperature for slave-boson-like solutions. 

FIG. 2. Phase Diagram of t/J=3 corrected by Gutzwiller's method. The labeling is the same 
as that in Fig.l. 
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